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Abstract 

By using the algebraic construction outlined in cni, we introduce several Markov 
processes related to the '^q(su(l,l)) quantum Lie algebra. These processes serve as 
asymmetric transport models and their algebraic structure easily allows to deduce duality 
properties of the systems. The results include: (a) the asymmetric version of the Inclusion 
Process, which is self-dual; (b) the diffusion limit of this process, which is a natural 
asymmetric analogue of the Brownian Energy Process and which turns out to have the 
symmetric Inclusion Process as a dual process; (c) the asymmetric analogue of the KMP 
Process, which also turns out to have a symmetric dual process. We give applications of 
the various duality relations by computing exponential moments of the current. 
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1 Introduction 


1.1 Motivations 

Exactly solvable stochastic systems out-of-equilibrium have received considerable attention 
in recent days [271 ng EH iH 0 [5]. Often in the analysis of these models duality (or self- 
duality) is a crucial ingredient by which the study of n-point correlations is reduced to 
the study of n dual particles. For instance, the exact current statistics in the case of the 
asymmetric exclusion process is obtained by solving the dual particle dynamics via Bethe 
ansatz IMllISlIl]. 

The duality property has algebraic roots, as was first noticed by Schiitz and Sandow for 
symmetric exclusion processes [25], which is related to the classical Lie algebra su(2). Next 
this symmetry approach was extended by Schiitz [20] to the quantum Lie algebra '^q(su(2)) in 
a representattion of spin 1/2, thus providing self-duality of the asymmetric exclusion process. 
Recently Markov processes with the '^q(su(2)) algebraic structure for higher spin value have 
been introduced and studied in [TOj. This lead to a family of non-integrable asymmetric 
generalization of the partial exclusion process (see also |22jL 

In [ig El] the algebraic approach to duality has been extended by connecting duality 
functions to the algebra of operators commuting with the generator of the process. In par¬ 
ticular for the models of heat conduction studied in m the underlying algebraic structure 
turned out to be '^(su(l, 1)). This class is richer than its fermionic counterpart related to 
the classical Lie algebra '^(su(2)) which is at the root of processes of exclusion type. In 
particular, the classical Lie algebra (su(l, 1)) has been shown to be related to a large class 
of symmetric processes, including: (a) an interacting particle system with attractive inter¬ 
actions (inclusion process [HI ES]); (b) interacting diffusion processes for heat conduction 
(Brownian energy process [Ill[9]); (c) redistribution models of KMP-type [T9[[8]. The dua¬ 
lities and self-dualities of all these processes arise naturally from the symmetries which are 
built in the construction. 

It is the aim of this paper to provide the asymmetric version of these models with (self)- 
duality property, via the study of the deformed quantum Lie algebra ^g(su(l,l)). This 
provides a new class of bulk-driven non-equilibrium systems with duality, which includes in 
particular an asymmetric version of the KMP model m- The diversity of models related to 
the classical '^(su(l, 1)) will also appear here in the asymmetric context where we consider 
the quantum Lie algebra '^g(su(l, I)). 

1.2 Models and abbreviations 

For the sake of simplicity, we will use the following acronyms in order to describe the class 
of new processes that arise from our construction. 

(a) Discrete representations will provide interacting particle systems in the class of Inclusion 
Processes. For a parameter k G M+, the Symmetric Inclusion Process version is denoted 
by SIP(/c), and ASlP{q,k) is the corresponding asymmetric version, with asymmetry 
parameter q G (0,1). 

(b) Continuous representations give rise to diffusion processes in the class of Brownian 
Energy Processes. For k G M_|_, the Symmetric Brownian Energy Process is denoted by 
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BEP(/c), and ABEP (cr, k) is the asymmetric version with asymmetry parameter a > 0. 

(c) By instantaneous thermalization, redistribution models are obtained, where energy or 
particles are redistributed at Poisson event times. This class includes the thermal- 
ized version of ABEP(cr, fe), which is denoted by Th-ABEP(cj,/c). In the particular 
case k = 1/2 the Th-ABEP((T,/c) is called the asymmetric KMP (Kipnis-Marchioro- 
Presutti) model, denoted by AKMP(cj), which becomes the KMP model as cr —)• 0. 
The instantaneous thermalization of the ASIP(g',/c) yields the Th-ASIP(g', A;) process. 

1.3 Markov processes with algebraic structure 

In [lOj we constructed a generalization of the asymmetric exclusion process, allowing 2j 
particles per site with self-duality properties reminiscent of the self-duality of the standard 
ASEP found initially by Schiitz in [26]. This construction followed a general scheme where 
one starts from the Casimir operator C of the quantum Lie algebra '^q(su(2)), and applies 
a coproduct to obtain an Hamiltonian working on the occupation number variables 

at sites i and i + 1. The operator H = then naturally allows a rich class 

of commuting operators (symmetries), obtained from the n-fold coproduct applied to any 
generator of the algebra. This operator H is not yet the generator of a Markov process. But 
H allows a strictly positive ground state, which can also be constructed from the symmetries 
applied to a trivial ground state. Via a ground state transformation, H can then be turned 
into a Markov generator L of a jump process where particles hop between nearest neighbor 
sites and at most 2j particles per site are allowed. The symmetries of H directly translate 
into the symmetries of L, which in turn directly translate into self-duality functions. 

This construction is in principle applicable to every quantum Lie algebra with a non¬ 
trivial center. However, it is not guaranteed that a Markov generator can be obtained. This 
depends on the chosen representation of the generators of the algebra, and the choice of the 
co-product. Recently the construction has been applied to algebras with higher rank, such 
as '^q(gl(3)) jS] [20] or '^g(sp(4)) [20], yielding two-component asymmetric exclusion process 
with multiple conserved species of particles. 

1.4 Informal description of main results 

In |14] we introduced a class of processes with su(l, 1) symmetry which in fact arise from this 
construction for the Lie algebra '^(su(l, 1)). In this paper we look for natural asymmetric 
versions of the processes constructed in [14] . and |8]. In particular the natural asymmetric 
analogue of the KMP process is a target. The main results are the following 

(a) Self-duality of ASIP{q,k). We proceed via the same construction as in [10] for the 
algebra '^g(su(l, 1)) to find the ASIP(g, k) which is the “correct” asymmetric analogue 
of the SIP(/i;). The parameter q tunes the asymmetry: q ^ 1 gives back the SIP(A:). 
This process is then via its construction self-dual with a non-local self-duality function. 

(b) Duality between ABEP{a, k) and SIP{k). We then show that in the limit e — >• 0 where 

simultaneously the asymmetry is going to zero {q = 1 — ea tends to unity), and the 
number of particles to infinity rji = we obtain a diffusion process ABEP(ct, A:) 

which is reminiscent of the Wright-Pisher diffusion with mutation and a selective drift. 
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As a consequence of self-duality of ASIP(g, k) we show that this diffusion process is dual 
to the SIP (A:), i.e., the dual process is symmetric, and the asymmetry is in the duality 
function. Notice that this is the hrst example of duality between a truly asymmetric 
system (i.e. bulk-driven) and a symmetric system (with zero current). 

(c) Duality of instantaneous thermalization models. Finally, we then consider instantaneous 
thermalization of ABEP((T,/c) to obtain an asymmetric energy redistribution model of 
KMP type. Its dual is the instantaneous thermalization of the SIP(A:) which for k = If 2 
is exactly the dual KMP process. 

1.5 Organization of the paper 

The rest of our paper is organized as follows. In section[2]we introduce the process ASIP(g, k). 
After discussing some limiting cases, we show that this process has reversible prohle product 
measures on Z_|_ (but not on Z). 

In section [3] we consider the weak asymmetry limit of ASIP(g',A;). This leads to the 
diffusion process ABEP(it, A:), that also has reversible inhomogeneous product measures on 
the half-line. We prove that ABEP(it, A:) is a genuine non-equilibrium asymmetric system in 
the sense that it has a non-zero average current. Nevertheless in the last part of section [3] 
we show that the ABEP((T, k) can be mapped - via a global change of coordinates - to the 
BEP(A:), which is a symmetric system with zero-current. In section IHlUl this is also explained 
in the framework of the representation theory of the classical Lie algebra ‘^(5u(l, 1)). 

In section m we introduce the instantaneous thermalization limits of both ASIP(q', k) and 
ABEP((T, j) which are a particle, resp. energy, redistribution model at Poisson event times. 
This provides asymmetric redistribution models of KMP type. 

In section[5]we introduce the self-duality of the ASIP(( 7 , k) and prove various other duality 
relations that follow from it. In particular, once the self-duality of ASIP(q,A:) is obtained, 
duality of ABEP((T, A;) with SIP (A;) follows from a limiting procedure which is proved in 
Section 15.21 In the limit of an inhnite number of particles with weak-asymmetry, the original 
process scales to ABEP((T, k), whereas in the dual process the asymmetry disappears because 
the number of particles is finite. Next the self-duality and duality of thermalized models is 
derived in Section 15.31 

In section [6] we illustrate the use of the duality relations in various computations of 
exponential moments of currents. Finally, the last section is devoted to the full construction 
of the ASIP(( 7 , k) from a ‘^q(su(l, 1)) symmetric quantum Hamiltonian and the proof of self¬ 
duality from the symmetries of this Hamiltonian. 

2 The Asymmetric Inclusion Process ASIP(g,/c) 

2.1 Basic notation 

We will consider as underlying lattice the finite lattice A/, = {1,... , L} or the periodic lattice 
Tl = T^jUL. At the sites of A^, we allow an arbitrary number of particles. The particle 
system conhguration space is Dl = Elements of are denoted by T],f^ and for ?? E 0^, 
i E Al, we denote by E N the number of particles at site i. For r] E and i,j E such 
that r]i > 0, we denote by the configuration obtained from rj by removing one particle 
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from i and putting it at j. 


We need some further notation of g-numbers. For q G (0,1) and n G Nq we introduce the 
g-number 



( 2 . 1 ) 


satisfying the property limg^i[n]g = n. The hrst g-number’s are thus given by 


[0]g = 0, [l]g = l, [2]g = q + q-\ [3], = q^ + 1 + q-\ 


We also introduce the g'-factorial 


[n]gl := [n]q ■ [n - l]g . [l]g , 


and the g-binomial coefficient 



Further we denote 


(a; q)m ■= (1 - a)(l - ag) • • • (1 - ag”" . 


( 2 . 2 ) 


m • — 


2.2 The ASIP(g,/c) process 

We introduce the process in finite volume by specifying its generator. 

Definition 2.1 (ASIP(q,k) process). 

1. The ASIP{q, k) with closed boundary conditions is defined as the Markov process onQ^ 
with generator defined on functions / : ^ M 


L-l 








+ + r?4[r?,+i],(/(V+i’‘) - /(r?)) (2.3) 


2. The ASIP{q, k) with periodic boundary conditions is defined as the Markov process on 
with generator 



(2.4) 


Since in finite volume we always start with finitely many particles, and the total particle num¬ 


ber is conserved, the process is automatically well defined as a finite state space continuous 
time Markov chain. Later on (see Section [6.1 II we will consider expectations of the self-duality 
functions in the infinite volume limit. In this way we can deal with relevant infinite volume 
expectations without having to solve the full existence problem of the ASIP(g, k) in infinite 
volume for a generic initial data. This might actually be an hard problem due to the lack of 
monotonicity. 
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2.3 Limiting cases 

The ASIP (q, k) degenerates to well known interacting particle systems when its parameters 
take the limiting values g —)• 1 and k ^ oo recovering the cases of symmetric evolution or 
totally asymmetric zero range interaction. Notice in particular that these two limits do not 
commute. 

• Convergence to symmetric processes 

i) q— )• l,k fixed: The ASIP(g',/c) reduces to the SIP(/c), i.e. the Symmetric Inclu¬ 
sion Process with parameter k. All the results of the present paper apply also to 
this symmetric case. In particular, in the limit q ^ 1, the self-duality functions 
that will be given in theorem 15.11 below converge to the self-duality functions of 
the SIP(/c) (given in [8]). 

ii) q —)• 1, k —)> oo: Furthermore, when the symmetric inclusion process is time changed 
so that time is scaled down by a factor 1/2A;, then in the limit k ^ oo the symmet¬ 
ric inclusion converges weakly in path space to a system of symmetric independent 
random walkers (moving at rate 1). 

• Convergence to totally asymmetric processes 

hi) k —>• oo,q fixed: If the limit A: ^ oo is performed first, then a totally asymmet¬ 
ric system is obtained under proper time rescaling. Indeed, by multiplying the 
ASIP(g', k) generator by (1 — q‘^)q^^~^ one has 

(1 - in) = [/(„*’'+■) - /(,))] 

Therefore, considering the family of processes := labeled by 

k > 0 and defining 

one finds that in the limit k ^ oo the process y^^\t) converges weakly to the 
Totally Asymmetric Zero Range process y{t) with generator given by: 

^ ^ „—2yj_|_i 1 

= E „-2 _ 1 - /(y)]> ( 2 . 5 ) 

i=i ^ 

In this system, particles jump to the left only with rates that are monotone in¬ 
creasing functions of the occupation variable of the departure site. Note that the 
rates are unbounded for —>■ oo, nevertheless the process is well defined even 

in the inhnite volume, as it belongs to the class considered in [2]. This is to be 
compared to the case of the deformed algebra t/g(s[ 2 ) [TO] whose scaling limit with 
infinite spin is given by H) 

L—1 - 2y 

= E ^ ^ Sit ^ » (2-6) 

i=l ^ 
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Here particles jump to the right only with rates that are also a monotonous in¬ 
creasing function of the occupation variable of the departure site, however now it 
is a bounded function approaching 1 in the limit oo. In |12] it is proved that 
the totally asymmetric zero range process (12.6p is in the KPZ universality class. 
It is an interesting open problem to prove or disprove that the same conclusion 
holds true for (1^ [23]. We remark that the rates of (USD are (discrete) convex 
function and this also translates into convexity of the stationary current j{p) as a 
function of the density p, whereas for (|2.6I) we have concave relations. 

iv) k —>• oo, q —>• 1: In the limit q ^ 1 the zero range process in (j2.5l) reduces to a 
system of totally asymmetric independent walkers. This is to be compared to item 
ii) where symmetric walkers were found if the two limits were performed in the 
reversed order. 

2.4 Reversible profile product measures 

Here we describe the reversible measures of ASIP(g, k). 

Theorem 2.1 (Reversible measures of ASIP(g,/c)). For all L £ N, L > 2, the following 

results hold true: 

1.) the ASIP{q,k) on Al with closed boundary conditions admits a family labeled by a of 
reversible product measures with marginals given by 



(2.7) 


for i £ Al and a £ [0,g (with the convention (^^^ = Ij- The normalization 

is 



and for this measure 



2k-l 


(2.9) 


2.) The ASIP{q, k) process on the torus with periodic boundary condition does not admit 


homogeneous product measures. 


Proof. The proof of item 2.) is similar to the proof of Theorem 3.1, item d) in [ID] and we 
refer the reader to that paper for all details. To prove item 1.) consider the detailed balance 
relation 



where the hopping rates are given by 


+ ry+i]. 
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Cg(r?*’*+\r?) = + ry, - l],[ry,+i + 1], 

and ^ denotes a reversible measure. Suppose now that /r is a product measure of the form 
/i = Then (|2.10p holds if and only if 

yLii(77j-l)//i+i(?7i+i + l)g“^^[2A:+ryi-l]q[ryi+i + l]q = yxj(?7i)//i+i(ryi+i)g^^[ryi]q[2A:+ryi+i]g (2.11) 


which implies that there exists a G M so that for all i G 


lJLi{n) 


Hi{n - 1 ) 


= aq 


4:ki 


[2k + n — 1]^ 


nl 


( 2 . 12 ) 


Then (j2.7p follows from (j2.12p after using an induction argument on n. The normalization 
is computed by using Corollary 10.2.2 of [T]. We have that 


< oo if and only if 0 < a < q i) ^ny i G A/, 


(2.13) 


As a consequence (since q < I and i = 1 is the worst case) a must belong to the interval 
[0, The expectation (|2.9I) is obtained by exploiting the identity 

□ 

The following comments are in order: 

i) vanishing asymmetry: in the limit g —)• 1 the reversible product measure of ASIP(g', A:) 
converges to a product of Negative Binomial distributions with shape parameter 2k and 
success probability a, which are the reversible measures of the SIP(A:) [8]. 

ii) monotonicity of the profile: the average occupation number E(“^(ryj) in formula (|2.9|) is 
a decreasing function of i, and limj_j.oo E(“)(77j) = 0. 

iii) infinite volume: the reversible product measures with marginal (1^ are also well- 
defined in the limit L —)• oo. One could go further to [—M, oo) H Z for a < gi4fcM+2fc-i 
(but not to the full line Z). These infinite volume measure concentrate on configurations 
with a finite number of particles, and thus are the analogue of the profile measures in 
the asymmetric exclusion process m- 


3 The Asymmetric Brownian Energy Process ABEP {(j,k) 

Here we will take the limit of weak asymmetry g = 1 — ecr —>■ 1 (e ^ 0) combined with the 
number of particles proportional to e~^, going to infinity, and work with rescaled particle 
numbers Xi = [er/jj. Reminiscent of scaling limits in population dynamics, this leads to a 
diffusion process of Wright-Fisher type [S], with cj-dependent drift term, playing the role 
of a selective drift in the population dynamics language, or bulk driving term in the non¬ 
equilibrium statistical physics language. 
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3.1 Definition 

We define the ABEP(g, A:) process via its generator. It has state space = (1^+)^) 1R+ := 
[0,+oo). Configurations are denoted by x G with Xi being interpreted as the energy at 
site i G hi. 


Definition 3.1 (ABEP(fT,/c) process). 

1. Let u > 0 and k > 0. The Markov process ABEP{a,k) on the state space ^i with 
closed boundary conditions is defined by the generator working on the core of smooth 
functions f : hTi —>■ M via 


L-l 




(3.1) 


2 = 1 


with 


^ (1 - e-2~,)(e2™.+. _ 1) 




d 


d 


dxi dxi+i 


fix) 


1 


|(1 - - 1) + 2A: (2 - fix) 

J \UXi UXi-\.l J 


2. The ABEP{a, k) with periodic boundary conditions is defined as the Markov process on 
with generator 


ieTz, 


The ABEP((T,/c) is a genuine asymmetric non-equilibrium system, in the sense that its 
translation-invariant stationary state may sustain a non-zero current. To see this, let E 
denote expectation with respect to the translation invariant measure for the ABEP((T, A;) on 
Ti. Let fiix) ■= Xi, then from (|3.2I) we have 

= 0i,j+i(x) - 0i_i,i(x) (3.3) 

with 

0i,i+i(x) = - ^ |(1 - e-2-^i)(e2-^2+i _ 1) + 2A) (2 - | (3.4) 

So we have 

[fiix{t))] = E,j; [0i,i+i(x(A))] - E,j; [0i_i,i(x(t))] 

and then, from the continuity equation we have that, in a translation invariant state, := 

—E [0i^i+i] is the instantaneous stationary current over the edge (i, i -|- 1). Thus we have the 
following 
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Proposition 3.1 (Non-zero current of ABEP(cr, A:)). 

= -E [0i,i+i] <0 ^/ A: > 1/2 

and 

= -E [0j,j+i] >0 if k = Q . 

Proof. In the case k > 1/2, taking expectation of (|3.4p we obtain 

E [0i,i+i] = ^ |(l - 4A:) + {2k - + ^{e^<^i^i+r-^i))^ 

Since expectation in the translation invariant stationary state of local variables are the same 
on each site and cosh(x) > 1 one obtains 


E 


[0i,+i]>^ {(1-4/^) + 2(2A;-1) + 


E 


^2a(xi+]_-Xi) 


Furthermore, Jensen inequality and translation invariance implies that 
E [0i,i+i] > ^ |(l - 4A:) + 2{2k - 1) + l} = 0 


In the case k = 0 one has 


1 


E[0,,,+i] = — E 

2a I 


< 0 


which is negative because the function is negative a.s. □ 


3.2 Limiting cases 

• Symmetric processes 


i) fj — )• 0, k fixed: we recover the Brownian Energy Process with parameter k, 
BEP(A:) (see [5]) whose generator is 







—) 

dxi+i J 


2 


2k{xi 



—i 

dxi+i J 


(3.5) 


ii) a —)• 0, k —7- oo: under the time rescaling t —?■ t/2k, one finds that in the limit 
k ^ oo the BEP(A:) process scales to a symmetric deterministic system evolving 
with generator 

( 3 , 6 ) 

This deterministic system is symmetric in the sense that if the initial condition is 
given by (xi(0), Xi+i(0)) = (a, b) then the asymptotic solution is given by the fixed 
point where the initial total energy a + 6 is equally shared among the 

two sites. 
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Wright-Fisher diffusion 


iii) cT ~ 0,k fixed: the ABEP(it, A:) on the simplex can be read as a Wright Fisher 
model with mutation and selection, however we have not been able to find in the 
literature the specific form of selection appearing in (|3.2p (see [S] for the analogous 
result when u = 0). To first order in a one recovers the standard Wright-Fisher 
model with constant mutation k and selection a 


^WF(a,k) 





—] 

dxi+i J 


2 


{2axiXi+i + 2k{xi 



—) 

dxi+i J 


• Asymmetric Deterministic System 


iv) k ^ oo, a fixed: if the limit A: —>■ oo is taken directly on the ABEP((T, k) then, by 
time rescaling t —?• t/2k one arrives at an asymmetric deterministic system with 
generator 


1 (o — _ f,2(7Xi+i\ f _ ^ 

This deterministic system is asymmetric in the sense that if the initial condition 
is given by (a:j(0), Xi+i(0)) = {a,b) then the asymptotic solution is given by the 
fixed point 



(A,B) 



1 _)_ g2(7(a-|-fe) 


, a + 6 - — In 
2i(7 


1 /1 + 


where A > B. 

v) k ^ oo, cj ^ 0: in the limit cr —>• 0 dazD converges to (ITGD and one recovers 
again the symmetric equi-distribution between the two sites of DEP process with 
generator dMl). 

vi) k ^ oo, cj —?■ oo: in the limit a ^ oo one has the totally asymmetric stationary 
solution (a + 6,0). 


3.3 The ABEP((j, A:) as a diffusion limit of ASIP(q', A:). 

Here we show that the ABEP(cj, k) arises from the ASIP(g, k) in a limit of vanishing asym¬ 
metry and infinite particle number. 

Theorem 3.1 (Weak asymmetry limit of ASIP(g, A;)). Fix T > 0. Let : 0 < T} 

denote the ASIP{1 — ae,k) starting from initial condition r}^{0). Assume that 

lim er7^(0) = x G (3.8) 

Then as e ^ Q, the proeess {r]’^{t) ■ 0 < t < T} eonverges weakly on path space to the 
ABEP{a, k) starting from x. 
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Proof. The proof follows the lines of the corresponding results in population dynamics 
literature, i.e., Taylor expansion of the generator and keeping the relevant orders. Indeed, 
by the Trotter-Kurtz theorem [21], we have to prove that on the core of the generator of the 
limiting process, we have convergence of generators. Because the generator is a sum of terms 
working on two variables, our theorem follows from the computational lemma below. □ 


Lemma 3.1. If rj^ G Hl is such that —>■ a: G then, for every smooth function F : 

—>■ IK; cmd for every i G {1,..., T — 1} we have 


where Ff:{rj) = F{erj), r] G Hl- 

Proof. Define = erj'^. Then we have, by the regularity assumptions on F that 


(3.9) 




= e 


d d 


dxi+i dxi 




d d 


dxi dxi+i 


F{x^) + 0{e^) (3.10) 


and similarly 




= —e 


d d 




d d 


dxi dxi+i 


dxi+i dxi 

Then using q = 1 — ea, and 

(1 - - 2xia^e-^^^^e + O(e^) 

straightforward computations give 


F(x^) + 0(e^) (3.11) 


■^2,2 


with 




B,(x^) 


d d 


dxij^i dxi 


+ D^{x^) 


d d 


dxi dxi+i 


F{x^) + 0{e) 


B,{x) = ^ |(1 - e-2-^*)(e2-^»+i - 1) + 2A; (2 - | + 0(e) 


D,{x) = ^ (1 - e-2--^)(e2--^+i - 1) + 0(e) 


(3.12) 


Then we recognize 




d d 


^dxi+i dxi 

= (x^) 


“ ^ + Beix'') - 


d d 


dxi dxi+i 


F{x^) 
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which ends the proof of the lemma by the smoothness of F and because by assumption, 

—>■ X. □ 

The weak asymmetry limit can also be performed on the q-TAZRP. This yields a totally 
asymmetric deterministic system as described in the following theorem. 

Theorem 3.2 (Weak asymmetry limit of g-TAZRP). Fix T > 0. Let {y^{t) : 0 < T} denote 
the q^-TAZRP, := 1 — ae, with generator (ESI) and initial condition y^{0). Assume that 

lim e/(0) = y e (3.13) 

£->■0 

Then as e —>■ 0, the process {y^(t) ■ 0 < t < T} converges weakly on path space to the Totally 
Asymmetric Deterministic Energy Process, TADEP{a) with generator 

(3.14) 

initialized from the configuration y. 

Proof. The proof is analogous to the proof of Theorem 13.11 □ 


3.4 Reversible measure of the ABEP (cr, k) 

Theorem 3.3 (ABEP((T,/c) reversible measures). For all L ^ 'H,L > 2, the ABEP{q, k) 
on with closed boundary conditions admits a family (labeled by 'j > —Aak) of reversible 
product measures with marginals given by 


;,.(xi) (1 - G M+ (3.15) 

X- 

for i G Al and 

= — Beta (2ki + —, 2k') (3.16) 

2cj V 2cj / 

Proof. The adjoint of the generator of the ABEP(cj, k) is given by 


with 




d 


1 


2a \dxi+i 



i=l 


- - 1) / j 

-1) + 2k [{I - - (e^'^^^+i 


(3.17) 
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Let /i be a product measure with ii{x) = then in order for /r to be a stationary 

measure it is sufficient to impose that the conditions 

- {(1 - - 1) + 2yt [(1 - - l)] } ^^{x) = 0 

are satisfied for any — 1}. This is true if and only if 


/i'(xi) 2k-e - 2/i; 

—7—T ~ - n - \- cr = -7-r + za —-— 

Hi{xi) 1-e Hi+i{xi+i) e 2 o-xi+i _ X 


(3.18) 


for any Xi,Xj+i G The conditions (I3.18P are verified if and only if the marginals /ij(x) 
are of the form (|3.15l) for some 7 G M, is a normalization constant, and the constraint 
7 > —Aak is imposed in order to assure the integrability of //(•) on 3t^L. Thus we have proved 
that the product measure with marginal (|3.15p are stationary. One can also verify that for 
any / : R 

^ABEP j 

which then implies that the measure is reversible. □ 


= i(^ 

k- 


ABEP 


Tipf) 


Remark 3.1. In the limit <7 —)• 0 the reversible product measure of ABEP{a,k) converges to 
a product of Gamma distributions with shape parameter 2k and scale parameter I/ 7 , which 
are the reversible homogeneous measures of the BEP{k) |^. In the case a ^ 0 the reversible 
product measure of ABEP{a, k) has a decreasing average profile (see Proposition \4- 1\) . 

3.5 Transforming the ABEP(cr, fc) to BEP(/i:) 

In this subsection we show that the ABEP(it, A:), which is an asymmetric process, can be 
mapped via a global change of coordinates to the BEP(k) process which is symmetric. Here 
we focus on the analytical aspects of such cr-dependent mapping. In Section 13.61 we will show 
that this map induces a conjugacy at the level of the underlying su(l, I) algebra. This implies 
that the ABEP(( 7 ,A:) generator has a classical (i.e. non deformed) su(l,l) symmetry. This 
is a remarkable because ABEP(g',A;) is a bulk-driven non-equilibrium process with non-zero 
average current (as it has been shown in Proposition 13.11) and yet is generator is an element 
of the classical 5u(l, 1) algebra. 

Definition 3.2 (Partial energy). ITe define the partial energy functions Ei : -A R+, 

z G {1,..., L -|- 1} 

L 

Ei{x) :='^^X£, for i ^ Al and El^i{x) = 0. (3.19) 

e=i 

We also define the total energy E : -A R+ as 

E{x) := Ei{x). 
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Definition 3.3 (The mapping g). We define the map g : 


^-2aEi+l{x) _ ^-2crEi{x) 

g{x) := {gi{x))i(zK^ with gfix) := - — - (3.20) 

Notice that g does not have full range, i.e. Indeed 

E{g{x)) = ^ (l - ^ (3.21) 

so that in particular C {x G '■ E{x) < l/2a}. Moreover is a bijection from 

to g[^L\- Indeed, for z G g{t%^L\ we have 


{g ^{z))i = 


1 Z]jr=i+1 ^3 1 


(3.22) 


Theorem 3.4 (Mapping from ABEP((T,/c) to BEP(A:)). Let X{t) = {Xi(t))iQA^ be the 
ABEP{a,k) process starting from X(0) = x, then the process Z{t) := {Zi{t))i^\j^ defined by 
the change of variable Z{t) := g{X(t)) is the BEP{k) with initial condition Z{0) = g{x). 

Proof. It is sufficient to prove that, for any / : M+ smooth, x G and g defined 

above 

[i?®+i/] {g{x)) = o g)]{x) (3.23) 

for any i G Al- Define F ■= f o g, then 


[^^^^^ifog)]{x) = [^^“^^(F)](x) = 


ABEP/ 


(3.24) 


_ e-2axi^^^2axi+i _ 


4cj2 

1 


d 


dx. 


i+l 


A. 

dxi 


F{x) 


+ I (1 - - 1) + 2A: (2 - \ (—^- F{x) 

ZCT J \^UXi-\.l UXi J 


The computation of the Jacobian of g 


% 

dxj 


—2agj{x) for J <i — 1 
(x) = <( for j = i 

0 for J > f + 1 


(3.25) 


implies that 


and 


d 


dxj 


i+l 


A 

dx- 


gj{x) = < 


( 0 for j < i — 1 

— g 2(TEi^_i(a;) ^qy j — i 

^-2aEi+i{x) for J = i + 1 

0 for j > i + 2 


^ ) F{x) = 

dxi+i oxi 


—- — 1 / 

dzi+i dzi 


ig{x)) 


(3.26) 


(3.27) 
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(3.28) 


axi+i oxi 


_|_ g-4o-£;i+i(x) 


" -^ 1 / 


dzi+i dzi 


d d 


{9{x)) 


dzi+i dzi 


{9{x)). 


Then, using (j3.27l) and (|3.28p . (|3.24p can be rewritten as 

[^“r(/»9)iw = 

= ^ n _ p-2iTXi\/ 2(7Xi+i _ -4aEi+i(x) f _ ^ \ r 

4(72^ 


(5(a;)) 


+ <j 2^7 + ^ { (1 - - 


1 ) + 2 A: (2 - ^ W 


5 9 


9zi+i 


i9{x)) 


Simplifying, this gives 

if ° 9)]{x) = 


3-2(TEi+i(x) _ „-2aEi{x) „-2cri?i+2(x) _ p-2aEi+i{x) 


2a 2(7 

-2aEi{x) _ 2g-2o-£;i+i(x) _|_ g-2o-Ei+2(x) 

= [-^li+T/] ( 5 ( 2 ;)) 


a a 


fc 

-I e 

a 


dzi+i dzi 

d d 


{9{x)) 


dzi+i dzi 


f 


{9{x)) 


□ 


The ABEP(c 7 , k) has a single conservation law given by the total energy E{x) = 

As a consequence there exists an infinite family of invariant measures which is hereafter 
described. 


Proposition 3.2 (Microcanonical measure of ABEP(( 7 ,/c)). The stationary measure of the 
ABEP{a, k) process on with given total energy E is unique and is given by the inhomoge¬ 
neous product measure with marginals (|3.15l) conditioned to a total energy E{x) = E. More 
explicitly 


dpi^^\y) 


/■■■I nf=i y^=E}dyi 


(3.29) 


Proof. We start by observing that the stationary measure of the BEP(/c) process on A/, 
with given total energy S is unique and is given by a product of i.i.d. Gamma random variable 
(Xj)jgAA with shape parameter 2k conditioned to YlieAi ^■ This is a consequence of 

duality between BEP(/c) and SIP(/c) processes [13]. Eurthermore, an explicit computation 
shows that the reversible measure of ABEP(c 7 ,/c) conditioned to energy E are transformed 
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by the mapping g (see Definition ESD to the stationary measure of the BEP(A:) with energy 
S' given by 

The uniqueness for ABEP(it, k) follows from the uniqueness for BEP((T, k) and the fact that 
5 is a bijection from to g[^L\- D 


3.6 The algebraic structure of ABEP {a, k) 

First we recall from m that the BEP(A:) generator can be written in the form 


_^BEPik) ^ {K+Kr^, + K-Kt,, - + 2k^) (3.30) 

i=l 

where 

= 

K- = 

Kt = 

is a representation of the classical su(l,l) algebra. We show here that the ABEP((T,/c) has 
the same algebraic structure. This is proved by using a representation of su(l, 1) that is 
conjugated to ()3.31l) and is given by 


(3.31) 


^2 d 

o 2 2A: 

OZ^ OZi 

d 

^*■5 -^ ^ 

UZi 


kf = CgoKfo Cg-1 with a G {+, o} (3.32) 

where g is the function of Definition 13.31 and 


{Cg-if){x) = {f o g l)(x) 

iCgf){x) = if og){x) . 

Explicitly one has 

{Kif)ix) = (Af / o g-^){g{x)) with a € {+, o} (3.33) 


Theorem 3.5 (Algebraic structure of ABEP((T,/c)). The generator of the ABEP{a,k) pro¬ 
cess is written as 

L-l 

^ABEP[a,k) ^ ^ + k-kl, - kfkt^, + 2k^) (3.34) 

i=l 


where the operators kf with a G {+,—, 0 } are defined in (|3.32l) and provide a representation 
of the 5u(l, 1) Lie algebra. 
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Proof. The proof is a consequence of the following two results: 

^ABEP{a,k) _^BEP{k) ^ ^3 35 ^ 

and the operators Kf' with a E {+,—, 0 } satisfy the commutation relations of the 5u(l,l) 
algebra. The first property is an immediate consequence of Theorem 13.41 as Eq. (|3.35p is 
simply a rewriting of Eq. (I3.23D by using the definition of Cg and Cg-i. The second property 
can be obtained by the following elementary Lemma, which implies that the commutation 
relations of the iL“ operators with a E {+,—, 0 } are the same of the iL? operators with 
a E {+, —, o}. D 

Lemma 3.2. Consider an operator A working on function f : 1^ and let g : 

X C be a bijection. Then defining 

A = Cg 0^0 Cg—1 

we have that A^ A is an algebra homomorphism. 

Proof. We need to verify that 

A + B = A + B and AB = AB 

The first is trivial, the second is proved as follows 

AB = CgoABo Cg-i = {CgoAo Cg-i) o{CgoBo Cg-i) = AB 


As a consequence 

[^] = [A,B]. 

□ 


4 The Asymmetric KMP process, AKMP(cr) 

4.1 Instantaneous Thermalizations 

The procedure of instantaneous thermalization has been introduced in [H]. We consider a 
generator of the form 

= (4.1) 

i 

where is such that, for any initial condition {xi,Xi+i), the corresponding process con¬ 

verges to a unique stationary distribution 

Definition 4.1 (Instantaneous thermalized process). The instantaneous thermalization of 
the process with generator AC in (SH) is defined to be the process with generator 

i 
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where 


■^1,1+1 f 


lim “ /) 

t—>-oo 


(4.2) 



1) yi) Vi+lj 3;j-|_2) ■ ■ ■ ) S^l) f • • • ) ^L^]d^(^xi,Xi^i) (.Ui^ Vi+l) 


In words, in the process with generator £/ each edge (i, i + 1) is updated at rate one, 
and after update its variables are replaced by a sample of the stationary distribution of the 
process with generator starting from (rcj, Xj+i). Notice that, by definition, if a measure 

is stationary for the process with generator then it is also stationary for the process 

with generator 


An example of thermalized processes is the Th-BEP(A:) process, where the local redistribution 
rule is 

{x,y) ^ {B{x + y),{l-B){x + y)) (4.3) 

with B a Beta(2/c, 2/;;) distributed random variable [9j. In particular for A; = 1/2 this gives 
the KMP process m that has a uniform redistribution rule on [0,1]. Among discrete models 
we mention the Th-SIP(A:) process where the redistribution rule is 

{n,m) ^ {R,n + m — R) (4.4) 

where R is Beta-Binomial(n + m,2k,2k). For k = lj2 this corresponds to discrete uniform 
distributions on {0,1,... , n + m}. Other examples are described in [9]. In the following we 
introduce the asymmetric version of these redistribution models. 


4.2 Thermalized Asymmetric Inclusion process Th-ASIP(g,/c) 

The instantaneous thermalization limit of the Asymmetric Inclusion process is obtained as 
follows. Imagine on each bond + to run the ASIP(( 7 , k) dynamics for an inhnite amount 
of time. Then the total number of particles on the bond will be redistributed according to 
the stationary measure on that bond, conditioned to conservation of the total number of par¬ 
ticles of the bond. We consider the independent random variables (Mi,..., Ml) distributed 
according to the stationary measure of the ASIP(g, k) at equilibrium. Thus Mj and Mj+i are 
distributed according to 


Pi^Hvi) ■= P^"^(Mi = rji) 


^ fr]i + 2k-l\ 

z<"'V 


and 


r/i G N (4.5) 


p!+U^*+i) •= = m+i) = 


frii+i + 2k-l 

(a) 

i+1 


z. 


hi+1 


■ q 


Ak{i+l)r]i+i 


Vi+l € N 
(4.6) 
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for some a € [0, (2A:+i)^_ Hence the distribution of Mj, given that the sum is fixed to 

Mi + = n + m has the following probability mass function: 


\ n + m) := P(Mj = r | Mj + Mj+i = n + m) 

+ m - r) 

^ N —Akr f 2/c — l\ 

= %,k{n + m)q J 

where r G N and ‘^q^kin + m) is a normalization constant. 


(4.7) 


2k + n + m — r — 1 
n + m — r 


Definition 4.2 (Th-ASIP((?,/c) process). The Th-ASIP{q,k) process on is defined as 
the thermalized discrete process with state space TLi and local redistribution rule 


{n,m) ^ {Rq,n + m — Rq) (4.8) 

where Rq has a q-deformed Beta-Binomial{n + m, 2k, 2k) distribution with mass function 
(BID. The generator of this process is given by 

L-1 Vi+Vi+l 

= ^ [fiVi,---,Vi-i,r,Vi + Vi+i-r,Vi+ 2 ,---,VL)-fiv)] + 

i=l r=0 

(4.9) 


4.3 Thermalized Asymmetric Brownian energy process Th-ABEP (a, k). 

We define the instantaneous thermalization limit of the Asymmetric Brownian Energy process 
as follows. On each bond we run the ABEP(it, k) for an infinite time. Then the energies on 
the bond will be redistributed according to the stationary measure on that bond, conditioned 
to the conservation of the total energy of the bond. If we take two independent random 
variables W and Aj+i with distributions as in (I3.15E i.e. 

f^fixi) := (1 - Xi G M+ (4.10) 

/ii+i(xi+i) := (1 - e-2^^«+i){2fc-i)e-(4^^h+i)+7)^»+i Xi+i G M+ (4.11) 

‘^i+l 

then the distribution of Xi, given the sum hxed to Xi + Xj+i = E, has density 


p{xi\Xi + Xi+i = E) = 


Pi{xi)p,i+i{E - Xj) 

ff fj.i(x)/ 2 i+i(E - x) dx 


= ^a,k{E) e^^*^** (1 - (l - e-2<"(®-^*)) 


1 2k-l 
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where ’^a,k{E) is a normalization constant. Equivalently, let Wi := XijE, then Wi is a 
random variable taking values on [0,1]. Conditioned to Xi + Xj+i = E, its density is given 
by 

c / \ ^ 2fc—1 

v.^k{w\E) = - l) (l - } (4.12) 

with 

p\ , / \ 2k—\ 

K,k{E) ■■= J - 1) (l - } dw (4.13) 

Definition 4.3 (Thermalized ABEP((T,/c)). The Th-ABEP{a, k) process on Kl is defined 
as the thermalized process with state space and local redistribution rule 


(x, y) {B„{x + y), (1 - B„){x + y)) 


(4.14) 


where B^ has a distribution with density function r'(j,fc(-|ic + y) in (I4.12p . Thus the generator 
of Th-ABEP{a, k) is given by 


fix) = 

L-l „i 

= / [fixi, ■■■, wixi+ Xi+i), (1 

*=1 -^0 


(4.15) 


w){xi + Xi+i), ...,xl)- fix)] n^,kiw\xi +Xi+i) dw 


In the limit cr —>• 0, the conditional density z^q+;.(-|E) does not depend on E, and for any 
E > 0 we recover the Beta(2A:, 2k) distribution with density 


^o+,kM^') 


1 

Beta(2/i:, 2k) 


[u;(l 


lu)] 


2k-l 


(4.16) 


Then the generator coincides with the generator of the thermalized Brownian 

Energy process Th-BEP(A:) defined in equation (5.13) of [5]. 


The redistribution rule with the random variable B^ in Dehnition 14.31 is truly asymmetric, 
meaning that - on average - the energy is moved to the left. 

Proposition 4.1. Let B„ be the random variable on [0,1] distributed with density (|4.12p . 
then IE[Bo-] > As a consequence and 1 — are not equal in distribution and for 
(Xi ,... ,Xl) distributed according to the reversible product measure y of ABEPia, k) defined 
in (I3.15p . we have that the energy profile is decreasing, i.e. 

>E^[Ai+i], Vf e{1,...,L-1} . (4.17) 

Proof. Let X = iXi,X 2 ) be a two-dimensional random vector taking values in dis¬ 
tributed according to the microcanonical measure of ABEP(cj, k) with fixed total energy 
B > 0, defined in (13.2911 . Then, from Definition 14.31 

iXi,X 2 ) = iEB„,Eil-B„)) with B<,~p,,fc(-|B) (4.18) 
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Then, as already remarked in the proof of Proposition[321 Z := g{X) with g{-) as in Definition 
m\ is a two-dimensional random variable taking values in g[^ 2 ] C ^2 and distributed ac¬ 
cording to the microcanonical measure of BEP(A;) with hxed total energy S’ = ^(1 — 

It follows from (j4.3ll that 


g{X) = {SB,S{l- B)) with 5 ~ Beta(2fe, 2A:) 
Then, by (13.2211 we have 


{l-B^)E={g-\Z))2 = ^ln 


1 


2a [1- 2a{l - B)S 


and therefore 




(4.19) 


(4.20) 


(4.21) 


Put 2aE = 1 without loss of generality, for simplicity. Then to prove that E[i?a-] > 1/2 we 
have to prove that 

E(l + ln(l-B(l-e-i)))>i 
Defining a = 1 — e~^ we then have to prove that 


E(-ln(l-aB))<- 


It is useful to write 


ln(l-aB) = £ 




n=l 


n 


and remark that for a Beta{a, a) distributed B one has 


n—1 


E(B") = H 


r=0 


a + r 
2a+ r 


So we have to prove that 


00 r, n—1 


, , . V—> LL 1 —r Ot T . 

V’(a,«) := X] T n 


, n 2a -I- r 

n=l r=0 


First consider the limit a ^ 00 then we find 

Im (p{a,a) = -{1- e 


n=l 


-1^ 


1 e 


= — In ( —h 

2 2 


Next remark when a = 0 the B is distributed like ido + hSi which gives 


2 "u 2 '- 


E(-ln(l-aB)) =--ln(e ) = - 


(4.22) 


0.379 < 1/2 
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Now we prove that xjj is monotonically decreasing in a. To see this notice that 

d a + r _ -r 
da2a + r {2a + r)^ 


So the derivative 


OO 71—1 


da 


V^(a,a) 


71=1 r'=0 


a''' 

n 


71—1 

n 

, r=0,r^r' 


a + r 




2a + r I {2a + r)^ 


< 0 


Therefore o) is monotonically decreasing in a and '(/'(a, o) < ^. Thus the claim E[i7o-] > 
1/2 is proved. 

Now let X = {Xi^X2) be a two-dimensional r.v. distributed according to the profile 
measure + defined in (13.1511 with L = 2 and with abuse of notation let [.So-lE] = E 
Then we can write X = {E Bf^,E{l — Bu)) where now E is a random variable. We have 


E^[X2] = E^[E^[X2| E]] =E^[E^[E(1-B,)| E]] =E^[Ei/,,fc[(l-.B,)| E]] 

< E^ [E [Ba\ E]] = E^ [E^ [Xil E]] = E^ [Xi] (4.23) 

The proof can be easily generalized to the case L >2, yielding (I4.17p . □ 


For k = 1/2 and u —0 the Th-ABEP(cj,/c) is exactly the KMP process [19]. For A: = 1/2 
and a > 0 

^c7,i/2{w\E) =, u;G[0,1] (4.24) 

The Th-ABEP(cj, ^) can therefore be considered as the natural asymmetric analogue of the 
KMP process. This justifies the following definition. 

Definition 4.4 (AKMP((t) process). We define the Asymmetric KMP with asymmetry 
parameter a G M+ on Kl as the process with generator given by: 

^AKMP{a) + Xj+l) _ 

1=1 

[f{xi,w{xi + Xi+i), (1 - w){xi + Xi+i), ...,xl)- fix)] dw 

5 Duality relations 

In this section we derive various duality properties of the processes introduced in the previous 
sections. We start by recalling the definition of duality. 

Definition 5.1. Let {Xt]t>o, {Xt/tPO be two Markov processes with state spaces Ll and Ll 
and H : n X 0 —)• R a bounded measurable function. The processes {Xt}t>o, {Xt}t>o said 
to be dual with respect to D if 

E^[D{Xt,x)]=%[D{x,Xt)] (5.1) 
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for all X G Q,x & and t > 0. In a lEj; is the expectation with respect to the law of 
the {Xt}t>o process started at x, while Ej denotes expectation with respect to the law of the 
{Xt}t>o process initialized at x. 


5.1 Self-duality of ASIP(g, A;) 


The basic duality relation is the self-duality of ASIP(g, fc). This self-duality property is 
derived from a symmetry of the underlying Hamiltonian which is a sum of co-products of the 
Casimir operator. In [TO] this construction was achieved for the algebra '^g(su( 2 )), and from 
the Hamiltonian a Markov generator was constructed via a positive ground state. Here the 
construction and consequent symmetries is analogous, but for the algebra '^q(su(l, 1)). For 
the proof of the following Theorem we refer to Section 17.31 where we implement the steps of 
uni for the algebra '^q(su(l, 1 )). 

Theorem 5.1 (Self-duality of the finite ASIP(g, A;)). The ASIP{q,k) on Al with closed 
boundary conditions is self-dual with the following self-duality function 

L (Vi\ 










or, equivalently, 




2=1 




Jii-4:ki-e2Ni+i{r]))^i . 1 


with (a; q)m as defined in 


and 


(5.2) 


(5.3) 


^iiv) ■='^Vk 


k=i 


(5.4) 


Remark 5.1. For n G N, let ^/jg configurations with n particles located at sites 

Then for the configuration with one particle at site i 


„-(4fc£+l) 


(5.5) 


and, more generally, for the configuration ^(TvA) particles at sites with 

ii 7 ^ ij 

—4fcy^"_ t —IT? ^ 

Di,,,!'' .'■>)=■ n 


_ q—2k^n 


m=l 


The duality relation with duality function (|5.3p makes sense in the limit L —>• oo. Indeed, 
if Ni{rj) = oo for some i, then T)(^L'j{r], = 0 for all ^ with 7 ^ 0. If the initial 

configuration p G has a finite number of particles at the right of the origin, then from the 
duality relation, we deduce that it remains like this for all later times t > 0 , which implies 
that Nfirjt) < 00 for all t > 0. Conversely, if rj is such that No(rj) = 00 , then Ao(^t) = 00 for 
all later times because, from the duality relation, E^ [D^r], ^j)] = 0 for all t > 0. To extract 
some non-trivial informations from the duality relation in the infinite volume case, a suitable 
renormalization is required (see Section EH). 


25 






5.2 Duality between ABEP(cr, A:) and SlP(fc) 

We remind the reader that in the limit of zero asymmetry q ^ 1 the ASIP(g', k) converges to 
the SIP{k). Therefore from the self-duality of ASIP(g', k), and the fact that the ABEP((T, k) 
arises as a limit of ASIP(g, k) with g —>• 1, a duality between ABEP(cj, k) and SIP(A:) follows. 


Theorem 5.2 (Duality ABEP(fT, k) and SIP(A:)). The ABEP{a, k) on Kl with closed bound¬ 
ary conditions is dual to the SIP{k) on Al with closed boundary conditions, with the following 
self-duality function 






r(2A:) 

T(2k -|- 


^-2aEi+i{x) _ ^-2aEi{x) 




(5.6) 


with Ei{-) the partial energy function defined in Definition \S.S[ 

Proof. The duality function in (15.6p is related to the duality function between BEP(/c) and 
SIP(/c), D^^^{x,ri) (see e.g. Section 4.1 of |8]) by the following relation 

= Dfj ^) igix ), v ) (5.7) 

where g{-) is the map defined in (13.31) . Thus, omitting the subscript (L) in the following, 
from (|3.35p we have 




(rc) 




and this proves the Theorem. □ 


(5.8) 


Remark 5.2. In the limit os n —>■ 0 one recovers the duality Zl|’^^(-, •) between BEP{k) and 
SIP{k). However it is remarkable here that for finite a there is duality between a bulk driven 
asymmetric process, the ABEP{a,k), and an equilibrium symmetric process, the SIP{k). In¬ 
deed, the asymmetry is hidden in the duality function. This is somewhat reminiscent of the 
dualities between systems with reservoirs and absorbing systems m, where also the source 
of non-equilibrium, namely the different parameters of the reservoirs has been moved to the 
duality function. 


The following proposition explains how Zl^^^(x,|) arises as the limit of ASIP(( 7 , k) self-duality 
function for q = 1 — N~^a, N oo. 


Proposition 5.1. For any fixed L >2 we have 


lim 

N^OO 




(5.9) 
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where denotes the self-duality function of ASIP{q,k) defined in (|5.3p and 

denotes the duality function defined in (I5.6p . 

Proof. Let 

L 

N := \rj\ :='^r]i, q=l 

then 


i=l 


a 

N' 


X := N ^1], 


D 


ASIP(g,fc), _A [Vi]q[r]i-l]q---[Vi-Ci + ^q _ -L) [2 E 

(L) W,?; \A[2k + ^i-l]g[2k + C,-2],...[2k]g ^ 


(L) 

Now, for any m 


(5.10) 

(5.11) 


[r]i-m]i_- = [Nxi-mfi_j^ 


= ^ + 0(A^"^)] 


N 


a 


sinh((TXj) + 0(1) 


hence 


On the other hand 


Ci-i 


I N 

[Nxi - m]i_^ = ( — sinh(fJXj) + 0(1) 


m=0 




[2k-\-m]i_^ = 2k-\-m-\-0{N ^) thus [2k-\-m]i_^ 


m=0 


r(2fc + c,) 

r(2A:) 


(5.12) 

(5.13) 

+0(A^"^) (5.14) 


finally, let ffi^) := 2 ^rn + f.i and gfi^) := 2 Y!m=i im + f.i - ^kif,i we have 

Nxifiii) _, , / a\ 9 {i) 




= e-<^^i/i(0+o(Ar-i), and = (^1 - =1 + 0(1V-^) 

(5.15) 

then (|5.9p immediately follows. □ 


5.3 Duality for the instantaneous thermalizations 

In this section we will prove that the self-duality of ASIP(g,/c) and the duality between 
ABEP(cj, fc) and SIP(A:) imply duality properties also for the thermalized models. 

Proposition 5.2. If a process {r]{t) : t > 0} with generator ££ = Yld=i -^*,*+1 *5 dual to a 
process {^(t) : t > 0} with generator ^ = Yld=i -^*,*+1 duality function D{-, •) in such 
a way that for all i 

(r?) = K.+iO(7?,.)](0 

then, if the instantaneous thermalization processes of pt, resp. both exist, they are each 
other’s dual with the same duality function D{-, •). 
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Proof. Let ^ resp. be the generators of the instantaneous thermalization of rjt, resp. 
then, from (14.2p we know that 

£/ = Y' = lim - I) 

^ ^ t—>00 

and _ 

£/ = Y' M,i+i = lim - I) 

ieAh 

where I denotes identity and where the exponential is the semigroup generated by 

in the sense of the Hille Yosida theorem. Hence we immediately obtain that 


-/)!?(•, 0 (r/)= (e‘^-+i-/)D(r?,.) (0 




which proves the result. □ 


As a consequence of this Proposition we obtain duality between the thermalized ABEP(g, k) 
and the thermalized SlP(fe) as well as self-duality of the thermalized ASIP(g', k). 


Theorem 5.3. 

a) The Th-ASIP{q, k) with generator (j4.9p is self-dual with self-duality function given by 

(ISSD. 

b) The Th-ABEP{a, k) with generator (|4.15p is dual, with duality funetion (j5.6p to the 
Th-SIP{k) in Al whose generator is given by 

= (5.16) 

L—i 

= ^ • • • ’ ^*-1’ A + ^*+1 “ A ^i+2, fif)] I fi + ^j+i) 

i=l r=0 

where \ n + m) is the probability density of a Beta-Binomial distribution of pa¬ 

rameters (n + m, 2k, 2k). 


Remark 5.3. For k = 1/2 (I5.16p gives the KMP-dual, i.e., the asymmetric KMP has the 
same dual as the symmetric KMP, but of course with different a-dependent duality function 
given by 


j^AKMP{a) 

^(L) 




n 


1 


g _ g 2(TEj^{x'j 


ii 


(5.17) 


6 Applications to exponential moments of currents 

The definition of the ASIP(g, fc) process on the infinite lattice requires extra conditions on 
the initial data. Indeed, when the total number of particles is infinite, there is the possibility 
of the appearance of singularities, since a single site can accommodate an unbounded number 
of particles. By self-duality we can however make sense of expectations of duality functions 
in the infinite volume limit. This is the aim of the next section. 
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6.1 Infinite volume limit for ASIP(g,/c) 

In this section we approximate an infinite-volume configuration by a finite-volume configu¬ 
ration and we appropriately renormalize the self-duality function to avoid divergence in the 
thermodynamical limit. 


Definition 6.1 (Good infinite-volume configuration). 

a) We say that rj G is a “good infinite-volume configuration” for ASIP{q,k) iff for 
fj(L) g L G N, the restrietion of rj on {—L,L\, i.e. 


JL) ^ i Vi for i G [-L, L] 
* I 0 otherwise 


( 6 . 1 ) 


the limit 


lim 1 I g 
L—>-oo 




Ei: 




( 6 . 2 ) 




exists and is finite for allt > 0 and for any G finite (i.e. such that < ^)- 

b) Let be a probability measure on then we say that it is a “good infinite-volume 
measure” for ASIP{q, k) iff it concentrates on good infinite-volume configurations. 


Proposition 6.1. 


1) If ri G is a “good infinite-volume configuration” for ASIP{q,k) and 
configurations with n particles located at sites then the limit 


lim IT q 


L —^oo 


m=l 


^(^(i)^^(^i,...A)) 


(6.3) 


is well-defined for all t > 0 and is equal to 

n 


m=l 


(6.4) 


2) If 1 ] ^ is bounded, i.e. sup^g^??* < oo, then it is a “good infinite-volume configura¬ 
tion”. 

3) Let us denote by Plxit) a Poisson process of rate X > 0, and by E[-] the expectation w.r. 
to its probability law. If p, is a probability measure on such that for any A > 0 the 
expectation 

I Vi+i 


IE„ 


E 


(6.5) 


is finite for all t > 0 and for any i G Z, then fi is a “good infinite-volume measure’ 

Proof. 
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1) If r? G is a good infinite volume configuration, then the duality relation with duality 
function ()5.3p makes sense after the following renormalization: 






m=l 


m=l 


( 6 . 6 ) 

then the first statement of the Theorem follows after taking the limit as L ^ oo of 

dST]). 

2) Let ^ be a finite configuration in N^. We prove that for any bounded r] G the family 
of functions 


^L{t) := 


-2iiNi+-L(riW) 


Ef 




L G N 


(6.7) 


i£Z 


is uniformly bounded. Without loss of generality we can suppose that ^ ^ for 

some {£i, ... ,£n} C Z, n G N. Moreover we denote by ... ,in{'t)) the positions 

of the n ASIP(q,fe) walkers starting at time t = 0 from (^i,...,£„). We then have 
^(t) =^(LW,-A(t))^ and 


^ut) = n g 

m=l 




E 


= E, 






11 r ^ Ah . „9\ y 




. 2=1 


q-ikiru (t)+2[Ni^ (t)+i )-)] 




m=l 

As a consequence, since 
and 


1 


sup ■ 


for some c > 0, we have that there exists C > 0 such that 


U<n < 1 

(6.8) 

< c 

(6.9) 


|^L(t)| < C E^(,, 


jn=l 




( 6 . 10 ) 


for all L G N, t > 0. Then, from the Cauchy-Schwarz inequality, in order to hnd an 
upper bound for (|6.1Up . it is sufficient to hnd an upper bound for 


'SL,m(t) •— 


qK{-4:ki,nit)+2[N(^(t)+i{v^^^)-Ne^+l{ri(^'l)]} 

for any hxed m G {1,..., n} and k G N. Now, let M := supj^g rn < oo, then 
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hence there exists C',lj > 0 such that 


\sL,m{t)\ < C' 




( 6 . 11 ) 


for any L E N, t > 0. Since ^{t) has a finite number of particles, for each m E {1,..., n} 
the process \imit) — ^m\ is stochastically dominated by a Poisson process with 

parameter 

A:= max + 77 %} V max + r]]g[r]%} (6.12) 

0<7],rj' <n 0<r7,r?'<n 

then the right hand side of (| 6 . 1 ip is less or equal than 


E 




= e 


-At ^ Loi 


j=0 




(6.13) 


This proves that ^L{t) is uniformly bounded. 


3) Suppose that the probability measure /r satisfies (I6.5li . Then, in order to prove that it 
is a “good” measure, it is sufficient to show that 


lim Ep 

L—>-CxD 






.i£l. 


< 00 


(6.14) 


By exploiting the same arguments used in the proof of item 2), we claim that, in order 
to prove (16.1411 it is sufficient to show that for each fixed m = 1,... ,n, k > 0, the 
function 


■= Ep 


E 




K{-Akemit)+2[N,^ (T7(i) )-iVi^+i )]} 


(6.15) 


is uniformly bounded. We have that 


0 L,m(i) = 


= E„ 


< 


E 




E 








+ 1 


Then the result follows as in proof of item 2) from the fact that the process — is 
stochastically dominated by a Poisson process of rate A (|6.12jl . and from the hypothesis 

dSSD. 


□ 


Later on, if we write expectations in the infinite volume we always refer to the limiting 
procedure described above. Namely, for a “good infinite-volume configuration” 77 E N^, with 
an abuse of notation we will write 


TT q-^^im+iiv) := lim TT E . 7 .) 

L—^oo ' 


m=l 


m=l 


^(77(i),^(L,...A)) 

(6.16) 
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and 


n 5 [D{rj,^{t))] := ]J q Eg(«i,....«„) D{r]^^\^{t)) 

(6.17) 


m=l 


m=l 


6.2 ^-exponential moment of the current of ASlP{q, k) 

We start by defining the current for the ASIP(g, k) process on Z. 

Definition 6.2 (Current). Let {f?(t), t > 0} be a cadlag trajectory on the infinite-volume 
configuration space then the total integrated current Jfit) in the time interval [0,t] is 
defined as the net number of particles crossing the bond {i — l,i) in the right direction. 
Namely, let (ti)igN be the sequence of the process jump times. Then 

'7i(t) — ^ (6.18) 

k:tk£[0,t] 

Lemma 6.1 (Current). The total integrated current of a cadlag trajectory {r]{s))o<s<t with 
r/(0) = rj is given by 

Ji{t) = Ni{r]{t)) - Nfirj) := (^Ni{r]^^\t)) - (6.19) 

where Ni{g) is defined in (15.41) and is defined in (|6.1I) . Moreover 

lim Ji{t) = 0 (6.20) 


Proof. (I6.19p immediately follows from the definition of Ji{t), whereas (|6.20p follows from 
the conservation of the total number of particles. □ 


Proposition 6.2 (Current q-exponential moment via a dual walker). Let rj E a good 
infinite-volume configuration in the sense of Definition 1 6. il then the first q-exponential mo¬ 
ment of the current when the proeess is started from p at time t = 0 is given by 




qM(t) 


= q2iNiv)-Niiv)) 


i—1 

q-ikm{t) _ q-2vm(t)^ q‘i{Nm(t){v)-N{v)) 

n=—oo 




( 6 . 21 ) 

where m{t) denotes a continuous time asymmetric random walker on Z jumping left at rate 
q~‘^^[2k]q and jumping right at rate q‘^^[2k]q and Ej denotes the expectation with respect to 
the law of m{t) started at site i G Z at time t = 0. Furthermore N{r]) — Ni{r]) = 'Ylin<ihn o,nd 
the first term on the right hand side of (j6.2ip is zero when there are infinitely many particles 
to the left of i G Z in the eonfiguration rj. 


Proof. To prove (|6.21l) we consider the configuration E with a single dual particle 
at site i. Since the ASIP(( 7 , k) is self-dual the dynamics of the single dual particle is given an 
asymmetric random walk m{t) on Z whose rates are computed from the process definition 
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and coincides with those in the statement of the Proposition. Prom ()6.16|) . (|6.17p and item 
1) of Proposition 16.11 we have that 


is equal to 
q-^NM 

Then from (I6.19P we get 


-(Aki+l) 

<1 q-2NM ^ 


q2k _ q~2k 


'V 


q2Ni(v(t)) _ q2Ni+i{ri{t)) 


= q-2Ni{v) _i 


-1 


2k _ r,-2k 


q^r. _ q- 


q-ikmit) ^^2N^(t)(ri) _ q2N^q)+i{ri)'j 




^2Ji{t) 


= 


2Ji+i(t) 


+ Ej (^q^(^rn{t)('n)-Ni(ri)) _ ^(N^q)+i{ri)-Ni{ri))^^ (6.22) 


By iterating the relation in (I6.22p . for any n > 0 we get 


Ejy = g2(Afi-n(??)-Afi+l(f?)) E^ 


2Ji-„d) 


+ 


_'^q2{Ni_j{ri)-Ni+i{r]))^Ak{i-j) q-^km{t) ^^2{N^q){ri)-Ni_j{r])) _ ^2(7V^(j)+i (??)) ^ 

1=0 

(6.23) 


By taking the limit n —?> oo we get 


E^ 


2Ji+i{t) 


= lim q^i^i-nM-Ni+iiv)) E^ 


n—)-oo 
oo 


2Ji-n(t) 


+ 


_ y~^ g-2A^i+l(r?)g4fc(i-j) 

1=0 

and using (I6.20p we obtain ()6.2ip . □ 

We continue with a lemma that is useful in the following. 

Lemma 6.2. Let x{t) he the random walk on Z jumping to the right with rate a > 0 and to 
the left with rate b > 0, let a ^ M, and ACM then 


lim — log En 
t^oo t 


Oi 


,x{t) 


with 


x{t) G A = sup {x log a — J^{x)} — inf J^{x) 
J x£A 


+ y/x'^ + 4:ab 


^{x) = {a + h) — y/ x"^ + Aab + x In 


2a 


(6.24) 


(6.25) 
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Proof. Prom large deviations theory m we know that x{t)/t, conditioned to x{t)/t G A, 
satisfies a large deviation principle with rate function J^(x) — inix^A ^(x) where J^(x) is 
given by 

J^{x) := sup {zx — A(^;)} (6.26) 

Z 


with 


A(z) := lim - logE 

t—>co t 


,zx(t) 


= a(e^-l) + 6 (e-^-l) 


(6.27) 


from which it easily follows (j6.25D . The application of Varadhan’s lemma yields (|6.24p . 


□ 


Remark 6.1. Let m{t) he the random walk defined in Proposition \6.A then (|6.24D holds 
with 


J^{x) = [Ak]q - ^Jx"^ + {2[2k]qY + x log 


2[2k]qq 


2k 


X + ^ X'^ + {2[2k]qY 


(6.28) 


We denote by the expectation of the ASIP((;, k) process on Z initialized with the homo¬ 
geneous product measure on with marginals p, at time 0, i.e. 

E®^[/(?7(t))] = ^((8)i6z/i(r/i))E,,[/(r/(t))] . 

Proposition 6.3 (g-moment for product initial condition). Consider an homogeneous prod¬ 
uct probability measure fi on N. Then, for the infinite volume ASIP{q,k), we have 









q2Ji{t) 

— Eq 

V Ag j 

+ Eq 

q 4fcm(0 ^^rn{t) _ ^ 


(6.29) 

where \y := m{t) is the random walk defined in Proposition 16.SI In par¬ 

ticular we have 

lim - logE®^[g^'^*(*)] = sup {x log — J^{x)} — inf J^{x) (6.30) 

t^oo t x>0 x>0 


with Mq := q ^^\i/q and {x) given by (|6.28l) . 

Proof. It is easy to check that an homogeneous product measure p, verihes the condition 
(16.5p in Proposition 16.11 thus it is a good infinite-volume probability measure in the sense of 
Definition EH For this reason we can apply Proposition 16.21 and from (|6.2ip we have 






J ®p{dni) 


E^ 


2Ji{t) 


q 

i—1 


J (8)/r(dr/)g2(^W-^A^)) + ^ J (^q-2vm(t) _ i) q2(N„,g){n)-NM) 


Since 


34 
































then, in particular, f ^iiv)) = g since Xq < 1, where we recall the interpretation 

of N{r]) — Ni(r]) from Proposition 16.21 Hence 




with 


and 




i—1 „ 

{m{t) = m) 

i=-oo mez 


2{Nm+l{v)-NM) 


q 


— ~ l) M't) + (-^l/g - l) B{t) 

A{t) := (m(t) = m) 

n<i—l m<i 

B{t) := Y1 =m)q 


—4:km \ m—i 


1/9 


(6.33) 

(6.34) 


n<i—l m> 2 +l 


Now, let a := q then 


A{t) = ^ ^ P„ (m(t) = m) a 

n<i—l m<i 

j>l m<j 

= ^«-Po(^ 


m<0 


1-A, 




Xq Eq 


) = m) ^ Ai + ^ «”*Po {rn{t) = m 

i>i 

m>l 


^m(t) 1 

^ lm(t)<0 

+ Eq 

—4km(t) -t 

q I-m(t)>l 


j>m 


(6.35) 


Analogously one can prove that 

1 


B{t) = 


Ai/9 - 1 


|Eo /3”^Wi 


m(t)>2 


— Ai/qEo 


1^(,)>2] } (6.36) 


with /3 = g ^^Ai/g then (I6.29P follows by combining (16.321) . (16.351) and (I6.36|) . 

In order to prove ()6.30p we use the fact that m{t) has a Skellam distribution with parameters 
{[2k]qq‘^^t, [2k]qq~‘^^t), i.e. m{t) is the difference of two independent Poisson random variables 
with those parameters. This implies that 



/ \ 



Eq 

V A, j 

— Eq 

\m(t) 1 

'^q -^m(i)>0 


Then we can rewrite (j6.29p as 






— En 


\m(t) -I 


+ Po im{t) = 0) 


+ (A^ ^ — Ai/g) Eq q 


+ En 




— En 


lm(t)>0 (I + S’i{t) + S’2{t) + S’3{t) + S’4{t)) (6.37) 


-NiM) 

(6.32) 
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with 


^l{t) := 


Eq 


Eq 



Mt) ■= 


Po im{t) = 0) 


En 




and 


m) ■= 


Eq 

'\™(0i 


lm(t)>l] 

Eq 


Eq 




(6.38) 


To identify the leading term in (16.371) it remains to prove that, for each i = 1,2,3 there exists 
Cj > 0 such that 

sup < Cj (6.39) 

t>o 

This would imply, making use of Lemma 16.21 the result in (j6.30p . The bound in (|6.39p is 
immediate for i = 1,2, 3. To prove it for z = 4 it is sufficient to show that there exists c > 0 
such that 


A-'Eo 


—4km{t) -I 


(q 


m{t) 




A-^Eo 




< A-^Eo 


—4km(t) -I 


< cEo 

This follows since there exists m* > 1 such that for any m > m* A“^ < A™^ and then 

+ Eq 
+ Eq 

<(l + A-i)Eo' 

This concludes the proof. □ 


(6.40) 


< A-^Eo 






q-^^\ 


1/9 


m{t) 




(6.41) 


6.3 Infinite volume limit for ABEP {a, k) 

Definition 6.3 (Good infinite-volume configuration). 

a) We say that x E is a “good infinite-volume configuration” for ABEP{a,k) iff for 
E Ten, the restriction of x to [—L,L], i.e. 


xf) = 


Xi for i E [-L, L] 
0 otherwise 


the limit 


lim TT E. 

L^oo^ 


iez 


D%x^^\m) 


(6.42) 


(6.43) 


exists and is finite for aUt>0 and for any ^ E N^ finite (i.e. sueh that < ooj. 

b) Let jjL be a probability measure on then we say that it is a “good infinite-volume 
measure” for ABEP{a, k) iff it eoneentrates on good infinite-volume eonfigurations. 
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Proposition 6.4. 


1) If X & is a “good infinite-volume eonfiguration” for ABEP{a,k) and is the 

configurations with n particles located at sites ii,... ,in & then the limit 


m=l 

is well-defined for all t > 0 and is equal to 


lim 

L^OO 


n 


,2aEe^+l{xW) 




(L) 


lim IT 

^—>■00 ^ 
m=l 




(6.44) 


(6.45) 


2) If X ^ bounded, i.e. supjgg 

tion” for ABEP{a, k). 


Xi < oo, then it is a “good infinite-volume configura- 


3) Let us denote by Plxit) a Poisson process of rate A > 0, and by E[-] the expectation w.r. 
to its probability law. If p, is a probability measure on such that for any A > 0 the 
expectation 


E;, 


E 


(t) 

gZ^i=i 




(6.46) 


is finite for all t >0 and for any £ E Z, then p, is a “good infinite-volume measure” for 
ABEP{a,k). 


Proof. The proof is analogous to the proof of Proposition 16.11 □ 

Later on for a “good” infinite-volume configuration x E we will write 


H' 


,2(7^iEi+i{x) 


Eg[L»'"(x,C(t))] := lim 
L —>-00 


,2a^iEi+i{x(l^)) 


E, 




{X 


(L) 


,m) 


(6.47) 


and 

n 

L>'"(x(t),e^^i’-’^")) 

m=l 




n 



m=l 






L>'"(x(t),.e^^'’-’^"^) 

(6.48) 


6.4 e '^-exponential moment of the cnrrent of ABEP (a, k) 

We start by defining the current for the ABEP(it, k) process on Z. 

Definition 6.4 (Current). Let {x(t), t > 0} be a cddldg trajectory on the infinite-volume 
configuration space then the total integrated current Jfit) in the time interval [0, t] is 
defined as total energy crossing the bond {i — l,i) in the right direction. 

Ji{t) = Ei{x{t)) - Efixfif)) := lim (Efix^^^Xt)) - Ei{x^^'>)) (6.49) 

where Efix) is defined in (13.2h and x^^"> as in (I6.42p . 
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Lemma 6.3 (Current). We have limj^_(X) Jj(t) = 0. 

Proof. It immediately follows from the conservation of the total energy. □ 


Proposition 6.5 (Current exponential moment via a dual walker). The first exponential 
moment of Jfit) when the process is started from a “good infinite-volume initial configuration” 
X E at time t = is given by 








(6.50) 


where In{t) is the modified Bessel function. 


Proof. Let E be the configuration with a single particle at site £. Since the 
ABEP(cj,/c) is dual to the SIP(2/i;) the dynamics of the single dual particle is given by a 
continuous time symmetric random walker i{t) on Z jumping at rate 2k. Since x is a good 
configuration we have that the normalized expectation 


^2aEi{x) \D{x{t),C^^fi] = - E; 


Aka 


g-2crF£+l(a;(i)) _ ^-2uEi(x{t)) 


and, from the duality relation (15.51) this is also equal to: 
^2aEi{x) Z)(x,e(^(^))) 


o-2o-F£(t)_,.i(a;) _ -2(jEi(t){x) 


where denotes the expectation with respect to the law of i{t) started at site £ E Z at time 
t = 0. As a consequence, for any i G Z 


^2aEi(x) 


^-2aEi+i{x(t)) 


from which it follows 


^2aEi{x) jg 


^-2aEi{x{t)) 


= E, 


— g2o-Fi (x) 


^-2aEi(x{t)) 




=-2o-E^(t)+l(x) _ -2aEi(t){x) 


(6.51) 


E El 


t<i-l 


— g2o-Fi (x) 


E E. 

i<i-l 

= ^Eo 


,-2aEi(^t)+W) _ p-2(TEm)(x) 


^-2o-F£(i)+^+l(a:) _ -2aE^^^t)^t.{x) 


o-2crEi(^t) + m (^) 


E Eo 


m<i 


£< 2-1 


^-2aEi(^t)+d^) 


— fi2oEi(x) g-2o-F^(t)+i(x) 


= Ej 


^-2aEi(t){^) 


Thus we have arrived to 


E, 


^-2aJi{t) 


= Bi 


o-2<^{Et.{t)(x)-Ei{x)) 


(6.52) 


(6.53) 
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and the result ()6.50p follows since 


with 




Pi{l{t) = n) 


¥{l{t) = n I £(0) = i) 


where In{x) is the modified Bessel function. □ 


(6.54) 


Remark 6.2. Let i(t) be a continuous time symmetric random walk on Z jumping at rate 
2k, then ()6.24l) holds with 


J^{x) = Ak — -s/x'^ + {AkY + X log 


1 

Ak 


X 


+ + {Ak) 



(6.55) 


We denote by the expectation of the ABEP((T, fc) process on Z initialized with the 
omogeneous product measure on with marginals fi at time 0, i.e. 


E®^[/(x(t))] = J i^i^zhidxi)) Ea;[/(x(t))] 


(6.56) 


Proposition 6.6 (Exponential moment for product initial condition). Consider a proba¬ 
bility measure pL on M'*'. Then, for the infinite volume ABEP{a, k), we have 




g 2crJi(t) 


Po m = 0] + Eo 


(Af + Ai«) 




(6.57) 


where X± := f /j,(dy)e^^‘^^ and £{t) is the random walk defined in B.emark \6.Sl In particular 
we have 

lim - logE®^[e“^'^'^*^*^] = sup{xlogA+ — J^(x)} — inf J^(x) (6.58) 

t—^cx) t 

with ^{x) given by ()6.55p . 


Proof. It is easy to check that an homogeneous product measure /i verifies the condition 
(I6.46P in Proposition l6.ll thus it is a good infinite-volume probability measure for ABEP(ct, k) 
in the sense of Definition 16.31 Thus we can apply Proposition 16.51 in particular from (|6.53l) 
we have 






J (8)^(dx)E2 


p-2o-(®£(t) {x)-Ei(x)) 




J ®p,{dx)Pti 

(^(^) = ^) / ®pL{d: 


-2a{En{x)-Ei{x)) 


(6.59) 
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Since 


I = A!.-" l{n<^} + K~' hn>^} 


(6.60) 


it follows that 

g(g)/i 




Y, p* m = n) A!_-" + 

n<i n>i-\-l 


= Ei 

= Eq 
= Eq 


A 


i—l{t) - 


^e{t)<i 


e(t)-i 


,-m- 


A^(*)l 


L£(t)<0 


+ ^, 

_l_ 1 

+ A_|_ 1 

_1_ X^l^l 1 

“T A I 1 




£(t)>l 


- “T A+ ^£{t)>l 


where the last identity follows from the symmetry of i{t). Then (I6.57p is proved. 
In order to prove (|6.58p we rewrite ()6.57l) as 




— En 


A I 1 


£(t)>0 


with 


Si{t) := 


En 


A^O) ^ yl(t) 




(1 + + S’2{t)) 

Po ix{t) = 0) 


En 


Ai(‘)l 


e(t)>o 


m) ■= 


En 


x^(0i 

A 1 1 




where for i = 1, 2 there exists Cj > 0 such that 


SUp|(fi(t)| < Ci 
t>o 


(6.61) 


(6.62) 


(6.63) 


This and the result of Remark 16.21 conclude the proof of (|6.58p . □ 


7 Algebraic construction of ASIP(g,/c) and proof of the self¬ 
duality 

In this section we give the full proof of Theorem 15.11 It follows closely the lines of [TO] however 
the algebra and co-product are different. 

7.1 Algebraic structure and symmetries 
The quantum Lie algebra '^g(su(l, 1)) 

For q € (0,1) we consider the algebra with generators ,K ~, satisfying the commutation 
relations 


[K+,K-] = -[2K% [A°, A±] = ±A± , 


(7.1) 
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where [•, •] denotes the commutator, i.e. [A, B] = AB — BA, and 


[2KX : = 


g2K° _ q-2K^ 


q-q 


-1 


(7.2) 


This is the quantum Lie algebra '^g(su(l, 1)), that in the limit (? —)• 1 reduces to the Lie 
algebra 5u(l, 1). The Casimir element is 


C = [KX[K^ - 1], - K+K- 

A standard representation of the quantum Lie algebra '^g(su(l, 1)) is given by 


(7.3) 


K+\n) = y ^[r/ + 2fc]q[r/+ l]g |n + l) 

K \n) = XHq[v + ^k-l]q |n-l) (7.4) 

K^\n) = [r] + k) \n) . 


/c G N. Here the collection of column vectors |n), with n G N, denote the standard orthonor¬ 
mal basis with respect to the Euclidean scalar product, i.e. |n) = (0,... , 0,1,0,... , 0)^ with 
the element 1 in the position and with the symbol ^ denoting transposition. Here and in 
the following, with abuse of notation, we use the same symbol for a linear operator and the 
matrix associated to it in a given basis. In the representation (izai) the ladder operators 
and K~ are one the adjoint of the other, namely 

{K+y = K- (7.5) 


and the Casimir element is given by the diagonal matrix 


C|n) = [k]q[k - l]q\n) . 

We also observe that the '^g(su(l, 1)) commutation relations in (I7.ip can be rewritten as 
follows 


qKoj^+ = q K+q^o ( 7 , 5 ) 

qKoj^- ^ q-^K-q^° 

[K+,K-] = -[‘^K\ 


Co-product structure 

A co-product for the quantum Lie algebra ^g{5u{l, 1)) is defined as the map A : '^g(su(l, 1)) —>■ 

1)) (g) '^q(su(l, 1)) 

A{kX = q~^° + q^° ® , 

A(A:°) = ®1 + 1®K^ . (7.7) 

The co-product is an isomorphism for the quantum Lie algebra I7q(s[2), i-e. 

[A(iL+), A{K-)] = -[2A{KXq, [A(K°), A(7L±)] = ±A(A:±) . (7.8) 
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Moreover it can be easily checked that the co-product satisfies the co-associativity property 

(A (g) 1)A = (1 (g) A)A . (7.9) 

Since we are interested in extended systems we will work with the tensor product over copies 
of the ^g(5u(l, 1)) quantum algebra. We denote by K^,K~, Kf, with i G Z, the generators 
of the copy. Obviously algebra elements of different copies commute. As a consequence 
of (j7.9|) . one can define iteratively A"- : '^g(su(l, 1)) —'^g(su(l, i.e. higher power 

of A, as follows: for n = 1, from (IZ2D we have 

A(Af) = 

A(K») = 

for n > 2, 

A^{Kf) = A"- 
a”(a:°) = A"- 

n times 

The quantum Hamiltonian 

Starting from the quantum Lie algebra ^g{5u{l, 1)) and the co-product structure we would 
like to construct a linear operator (called “the quantum Hamiltonian” in the following and 
denoted by for a system of length L) with the following properties: 

1. it is '^q(su(l, 1)) symmetric, i.e. it admits non-trivial symmetries constructed from the 
generators of the quantum algebra; the non-trivial symmetries can then be used to 
construct self-duality functions; 

2. it can be associated to a continuos time Markov jump process, i.e. there exists a 
representation given by a matrix with non-negative out-of-diagonal elements (which 
can therefore be interpreted as the rates of an interacting particle systems) and with 
zero sum on each column. 


iLf (8) (g) 

+ (7.10) 

(Kf) (g) q~^°+i + (g) 

(77°) (g) 1 + 1 (g) ... ® 1 ®i7°+i . (7.11) 


A natural candidate for the quantum Hamiltonian operator is obtained by applying the co¬ 
product to the Casimir operator C in (17.3p . Using the co-product definition (17.7p . simple 
algebraic manipulations yield the following definition. 


Definition 7.1 (Quantum Hamiltonian). For every L G N, L > 2, we consider the operator 
defined by 

"(1) E ‘ = E ifiS '+=(^) • 

i=l i=l 

where the two-site Hamiltonian is the sum of 


C(L) 


(gr^/c — q 2k^^q2k 1 — q (2k 1) ^ 



L times 


(7.13) 


42 



and 


(7.14) 


, i,i+i — 

^(L) — ^ 


l0A{Ci) ® 1 ( 8 ) • - ( 8 ) 1 

(z—1) times (L—i—1) times 


and, from (17.3p and (17.71) . 

A{Ci) = A{Kt)A{K-) - A([i^O],)A([ifO _ 

Explicitely 


(7.15) 


- ^ + qq~‘^^^ ® - (g + g"^)} (7.16) 

Remark 7.1. By specializing (|7.16l) to the representation (17.411 we get 


A(CO = q^°\K+(^Kff,+K-®K+^,+ 




(7.17) 


(g" + g-")((?"-i + g-("-i)) f Ko _j,o 

g I — q X 


2{q-q-^y 

(gfc _g-fc)(^fc-l _g-(fc-l)) 


‘^{q-q~^y 


-K^ 

g-^i+l — q ^i+l 


[q«i + ® I -A?,. 


Remark 7.2. T/ie diagonal operator C(/^) m (I7.13p /las 6een added so that the ground state 
|0)(L) •= is a right eigenvector with eigenvalue zero, i.e. 77(l)| 0)(/,) = 0 as it is 

immediately seen using dial). 

Proposition 7.1. In the representation dzaD the operator 77(l) is self-adjoint. 

Proof. It is enough to consider the non-diagonal part of 77(l)- Using ()7.5p we have 

® Krf^q-^?+^ + q^?R- ® 77+ 

= R-q^° 0 77+1 + K+q^^ ® q~^^+^ 77^i 

= q^°+^Rr ® 77+1^-^?+!-^ + q^°-^R+ ® 77" ig-^°+i+^ 

where the last identity follows by using the commutation relations (j7.6p . This concludes the 
proof. □ 


Basic symmetries 

It is easy to construct symmetries for the operator by using the property that the 

co-product is an isomorphism for the '^g(su(1,1)) algebra. 
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Theorem 7.1 (Symmetries of Recalling (I7.11j) . we define the operators 


:= ^ (g) • • • (g) q^i-i q~^i+^ (g) ... (g q~^° , 

i=l 
L 

:= = y^ 'l (g) • y (g) 1 ®K^ (g) 1 (g) • y (g) 1 . (7.18) 

('^— 1 ) times times 

They are symmetries of the Hamiltonian (j7.12p . i.e. 

[H(L).-'ffL)l = [^'(l).-^(l)l=0. (7.19) 

Proof. We proceed by induction and prove only the result for (the case is similar). 
By construction '■= are symmetries of the two-site Hamiltonian 77^2)' Indeed 

this is an immediate consequence of the fact that the co-product defined in (|7.8I) conserves 
the commutation relations and the Casimir operator commutes with any other operator 
in the algebra : 

[77(2),77± ] = [A(Ci), A(ii:f)] = A([Ci,i7f]) = 0 . 

For the induction step assume now that it holds [H= 0. We have 




zyO 

^(L) 


gl)l + (7.20) 

The first term on the right hand side of ()7.20p can be seen to be zero using (17.lip with i = 1 
and n = L — 1: 

Distributing the commutator with the rule \A,BC] = B\A,C] + [A,B]C, the induction hy¬ 
pothesis and the fact that spins on different sites commute imply the claim. The second term 
on the right hand side of p7.20p is also seen to be zero by writing 


[h 


{L) 


.7^(1)! = l'> 


L-l,L 
(L) > 


^(L-2) 


^ g^(i-2) A(iF±_^)] = 0 


□ 


Remark 7.3. In the case q = l, the quantum Hamiltonian in Definition \ 7. 1 \ reduces to the 
(negative of the) well-known Heisenberg ferromagnetic quantum spin chain 

L-l 

H{l) = Y. + 2k^) , (7.21) 

i=l 

with spins Ki satisfying the su(l, 1) Lie algebra. The symmetries of this Hamiltonian are 
given by 

L L 

Hl) = E E Ed = E ■ 

2=1 2=1 
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7.2 Construction of ASlP{q, k) from the quantum Hamiltonian 

In order to construct a Markov process from the quantum Hamiltonian -H(l)) we make use of 
the following Theorem which has been proven in [lOj . 

Theorem 7.2 (Positive ground state transformation). Let A be a |n| x |n| matrix with 
non-negative off diagonal elements. Suppose there exists a column vector := 5 G with 
strictly positive entries and such that Ag = 0. Let us denote by G the diagonal matrix with 
entries G{x,x) = g{x) for x G Then we have the following 


a) The matrix 


= G-^AG 


with entries 

^{x, y) = x,yGnxn (7.22) 

9{x) 

is the generator of a Markov process {Xt : t > 0} taking values on 
h) S commutes with A if and only if G~^SG commutes with ^. 
c) If A = A*, where * denotes transposition, then the probability measure y on LI 


fx{x) 


{gjx)? 


is reversible for the process with generator ^. 


(7.23) 


Now we apply item a) of Theorem 17.21 with A = 77(l)- At this aim we need a non-trivial 
symmetry which yields a non-trivial ground state. Starting from the basic symmetries of 
II[L) described in Section 17771 and inspired by the analysis of the symmetric case (g = 1), it 
will be convenient to consider the exponential of those symmetries. 

The ^-exponential and its psendo-factorization 

Definition 7.2 (g-exponential). We define the q-analog of the exponential function as 

rpTl 

exp^(x) := ^ ^ (7.24) 

n>o 

where 

1 - o” 

'■= 1^7 ( 7 - 25 ) 

Remark 7.4. The q-numbers in (I7.25P are related to the q-numbers in (12.ip by the relation 
{n}q2 = [n]qq^~^. This implies {n}q2 \ = [n]q! and therefore 

n>0 
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Proposition 7.2 (Pseudo-factorization). 
such that for L E N and S' E M 

Let {gi,...,gL} and {ki,...,kL} be 

operators 

ki9i = rgih 

for ie Al . 

(7.27) 

Define 



L 

:= with ;= ki 

. ki for i > 1 and = 1, 

(7.28) 

then 

ex.pfig^^'>) = eyipfigi) ■ exp 

rik^^^92) .exp^(fc(^ 

(7.29) 

Moreover let 



L 

g(L) g. /i(*+i)^ with h(*) := k-^ ■ 

■ ■ ■ ■ kj} for i < L and = 1, 

(7.30) 

then 

expfig^^'i) = e^pfigi h^‘^'>) • 

• • • • ex.pfigL-1 h(-^)) • ex.pfigL) 

(7.31) 


See [To] for the proof. 

The exponential symmetry 5^^ 

In this Section we identify the symmetry that will be used in the construction of the process 
ASIP(g, fc). To have a symmetry that has quasi-product form over the sites we preliminary 


define more convenient generators of the '2^q(su(l, 1)) quantum Lie algebra. Let 

E-.= q^°K+, F:=K-q-^° and K := q^^° (7.32) 

From the commutation relations (EH) we deduce that {E, F, K) verify the relations 

KE = q'^EK and KF = q-'^FK [E, F] = - ^ . ( 7 . 33 ) 

q-q ^ 

Moreover, from Theorem rm the following co-products 

A{Ei) := A{q^°) ■ A{K+) = Ei <S> 1 + Ki <S> E2 (7.34) 

A(Fi) := A(iLf) • A{q-^°) =Fi® + 1 ® F 2 (7.35) 


are still symmetries of 77(2)- In general we can extend (|7.34p and (I7.35P to L sites, then we 
have that 

= A(^-i)(g^?).A(^-i)(iL+) 

= + q^^°+^2K+ + ... + 

= 771-1- K 1 E 2 + K 1 K 2 E 3 + ... + Ki ■ ... ■ Kl-iEl (7.36) 
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p{L) 

= + ... + Kl_^q-n-^-^Kl ^ 

= ■ ...-Kl^ + --- + Fl-i-KI^ +Fl (7.37) 

are symmetries of H. If we consider now the symmetry obtained by g-exponentiating 
then this operator will pseudo-factorize by Proposition 17.21 

Lemma 7.1. The operator 

5+) :=exp^.(ii;(^)) (7.38) 

is a symmetry of .f7{L) ■ Its matrix elements are given by 



r]i + 2k- 1\ 
Ci + 2k-l)g 


(7.39) 


Proof. From (|7.33|) we know that the operators Ei,Ki, copies of the operators dehned in 
(I7.32p . verify the conditions (j7.27p with r = q^. As a consequence, from (j7.36h . (I7.38p and 
Proposition 17.21 we have 

S’+) = expg2(F;(^)) 

= exp 2 (F;i) • exp 2 (A:iF; 2 ) • • • exp 2 (All • • • Kl-iEl) 


= expq2 {q^°K +) • expg2 {q^^°q^^Ki ^) • • • 

_ c+ c+ c+ 

— -^2 ■ ■ ■ 


2 „ 

expg 2 [q 


(7.40) 


where 5"^^ := expg 2 ^q‘^^m=i has been defined. Using (I7.26p . we find 




1 / 2 ^*“^ K° 


ei>o 






(7.41) 




A>0 


li 


li 


where in the last equality we used (|7.4F Thus we find 

5+)|6,...,a) = 5+^2+... 5+16,...,a) 


(7.42) 


6 + /6 + + 6 ~ 1 


y, V 


E n 

U,^2,...,U>0i=l \ 

. ^^4te+fc+i)+2£4EUii(«-+U^ 16+^1,... ,?L + 4) 


form which the matrix elements in (j7.39D are immediately found. □ 
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Construction of a positive ground state and the associated Markov process ASEP(g, j) 

By applying Theorem 17.21 we are now ready to identify the stochastic process related to the 
Hamiltonian in (17.121) . 

We start from the state |0) = |0,... ,0) which is obviously a trivial ground state of 
We then produce a positive ground state by acting with the symmetry in (I7.38p . Using 
(I7.42P we obtain 


|g) = S+|0,...,0)= n 


2k + li-l 


/ a 


U,^2,...,U>0*=1 V ^ ' 1 

Following the scheme in Theorem 17.21 we construct the operator defined by 

G{L)\m, ■,VL){m, • • .,vl\s^\o,.. . ,o) 


(7.43) 


In other words is represented by a diagonal matrix whose coefficients in the standard 
basis read 


(-71,, nL\G\L) Ifl,....?l) = n J^ (7.44) 

Note that G( 2 ,) is factorized over the sites, i.e. 

(ryi, ... ...,^l) = <S)f=i{vi\Gi\Ci) (7.45) 

As a consequence of item a) of Theorem 17.21 the operator conjugated to via G'^^y 

i.e. 

= G-i (7.46) 

is the generator of a Markov jump process r](t) = (r/i(t), • • •, rjiit)) describing particles jump¬ 
ing on the chain The state space of such a process is given hy and its elements are 
denoted by r/ = (r/i,..., t/^), where rn is interpreted as the number of particles at site i. The 
asymmetry is controlled by the parameter 0 < < 1. 

Proposition 7.3. The action of the Markov generator := is given by 

L-l 

(_j^U)j)(^) = ^(^i j+i/)(7/) with 
i=l 

= g"*-"*+i+(2fc-i)[^.]^[2A) + r/,+i],(/(r?G+i)_/(^)) 

+ + V^Url^+lWW^^n - fib)) (7-47) 

Proof. From Proposition 17.II we know that = 77 (l), hence we have that the operator 

^{L) ■= G[l)H{l)G7\ is the transposed of the generator defined by ()7.46l) . Then we 

have to verify that the transition rates to move from r] io f for the Markov process generated 
by (17.471) are equal to the elements (^|77(x,)|7/). 
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Since we already know that is a Markov generator, in order to prove the result it is 

sufficient to apply the similarity transformation given by the matrix G'(l) defined in (|7.44l) 
to the non-diagonal terms of (I7.16p . i.e. q^i . We show here the computation 

only for the first term, being the computation for the other term similar. 

We have 

{C^,C^+l\G^G,+, • 

= {C^\G,q^°K+Gr^\v^) ® {C^+l\G,+lK-^,q-^l^G-^,h+l) (7.48) 

Using (17.441) and (17.4p one has 

{^,\G,q^?KtG-^h) = + 1) (7.49) 

and 

{^i+,\G,+lK-^,q-^?^^G-^,\v^+l) = _ i) (7.50) 

Multiplying the last two expressions one has 

(r/*+'l%)|r/) = (7.51) 

that corresponds indeed to the rate to move from ry to 7 ^*+^’* in (|7.47|) . This concludes the 
proof. □ 


Remark 7.5. From item c) of Theorem \7.Sl we have that the product measure defined 
by 

h{L){v) = {v\G‘fL)\v) (7.52) 


is a reversible measure of .if . Notice that it corresponds to the reversible measure p(«) 
defined in (1221) with the choice a = 1. 


7.3 Self-Duality of ASIP(g, k) 

The following Proposition has been proven in [ID] and it will be key to the proof of ASIP(g, k 
self-duality. 


Proposition 7.4. Let A = A* be a matrix with non-negative off-diagonal elements, and 
g an eigenvector of A with eigenvalue zero, with strictly positive entries. Let ££ = G~^AG 
be the corresponding Markov generator. Let S be a symmetry of A, then G~^SG~^ is a 
self-duality function for the process with generator F£. 

We now use Proposition 17.41 and the exponential simmetry obtained in Section 17.21 to 
deduce a non-trivial duality function for the ASIP(g, A:) process. 


Proof of lOli in theorem I5.ll From Proposition 17.11 we know that is self-adjoint, 
then, using Proposition 17.41 with A = FL{p^, G = G(^L'j given by ()7.44p and S = 5)*^^ given by 
(I7.39P it follows that 


<^—1 c + 1 


(7.53) 
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is a self-duality function for the process generated by Its elements are computed as 

follows: 


{vis+m f ^ ') ■ 

_ ■TT / /f/A + 2A; — 1\ / /2A: -|- T/j — 1\ /2A: -|- — 1\ 

■fiV + V m )S ), 

■ X]m=i(?>Ti+fc)-l-Ci-l-fc-l-l] —(2A:i—fc+l)(T;i+,ri) . i 

n[5. + 2t-iyfe-6],! " 

Since both the original process and the dual process conserve the total number of particles 
it follows that in (|5.2p is also a duality function. □ 
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